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Abstract 

Usually the effects of isotropic inhomogeneities are not seriously taken into account in the de¬ 
termination of the cosmological parameters because of Copernican principle whose statement is 
that we do not live in the privileged domain in the universe. But Copernican principle has not 
been observationally confirmed yet in sufficient accuracy, and there is the possibility that there are 
non-negligible large-scale isotropic inhomogeneities in our universe. In this paper, we study the 
effects of the isotropic inhomogeneities on the determination of the cosmological parameters and 
show the probability that non-Copernican isotropic inhomogeneities mislead us into believing, for 
example, the phantom energy of the equation of state, p = wp with w < —1, even in case that 
w = — 1 is the true value. 
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I. INTRODUCTION 


Usually, we believe that large scale isotropic inhomogeneities whose symmetry center 
coincide with our location are so small that we do not need to take into account them in the 
determination of cosmological parameters. This belief is based on the so called Copernican 
principle, which is one of the basic working hypotheses in the physical cosmology. Copernican 
principle states that we do not live in the privileged domain in the universe. If this is true, 
the isotropy around our location means that our universe is isotropic at every point, or 
equivalently, homogeneous and isotropic, and hence the observed high isotropy of the cosmic 


microwave background radiation strong 
universe in global sense (see, e.g. Ref. 


;ly suggests the high homogeneity and isotropy of our 
1]) which, at the same time, imply the smallness of 
the large scale inhomogeneities isotropic in terms of our location. 

The universe model with large isotropic inhomogeneities have been studied in the context 
of the alternative scenario to explain the observed distance-redshift relation without intro¬ 
ducing the dark energy components within the framework of general relativity. Although 
there are several severe observational constraints for the alternative scenario to the dark 
2 - 


energy 


10], we should note that this universe model has not been completely excluded 


yet. Furthermore, it should be noted that even if there are dark energy components, the 


existence of isotropic inhomogeneities may significantly affects observational results 
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19]. 


Denoting the energy density and the pressure of the dark energy by p L \ and p^, its equation 
of state is given by 

Pd = wpd (1-1) 

with pa > 0 and w < —1/3. The special case of w = —1 corresponds to the cosmological 
constant. We would like to stress that if the observational data really implies that w is less 


than - 
of state 


, it is worthwhile to study the effects of isotropic inhomogeneities to the equation 
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19]; w < — 1 may cause the causality violation. This is the main purpose of the 


present paper. 

In this paper, we focus on the effects of isotropic inhomogeneities on the estimate of how 
large amount of the dark energy in the universe and its equation of state. We consider the 
universe model filled with dust and a positive cosmological constant and assume that there 
are large-scale isotropic inhomogeneities which can be treated by the linear perturbation in 
the Friedmann-Lemaitre-Robertson-Walker (FLRW) universe model. This universe model 
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has two arbitrary functions. In order to specify the two arbitrary functions, we need two 
conditions. One of these conditions is that we have considered only the growing mode 
of linear perturbation. To specify one more function, we solve the inverse problem which 
is the reconstruction of the universe model from observational data. We adopt distance- 
redshift relation as observational data. We solve the inverse problem with observational 
data which is the same as that of the universe filled with non-relativistic matter and dark 
energy with various w. To estimate the systematic error in determining the energy density 
and the equation of state of the dark energy due to isotropic inhomogeneities, we assume 
that inhomogeneities are described by the power spectrum of the random Gaussian statistics, 
which is consistent with the inflationary universe scenario. Then, we derive the probability 
that isotropic inhomogeneities obtained by solving the inverse problem are produced in the 
inflationary universe. 

The organization of this paper is as follows. In Sec. II, we derive the basic equations 
for inhomogeneous and isotropic universe model. In Sec. Ill, we derive the null geodesic 
equations in order to relate the distance-redshift relation with the universe model. In Sec. IV, 
we derive the basic equations to determine the inhomogeneous and isotropic universe model 
from a given distance-redshift relation. We explain the numerical procedure in Sec. V, and 
then we show the numerical results in Sec. VI. We derive the probability of occurrence of 
systematic errors, in Sec. VII. Finally, Sec. VIII is devoted to the summary and discussion. 

In this paper, we adopt the sign conventions of the metric and Riemann tensor of Ref.[l] 
and the geometrized unit in which the speed of light and Newton’s gravitational constant 
are one. 


II. INHOMOGENEOUS ISOTROPIC UNIVERSE MODEL 


As mentioned in Sec]H we consider the FLRW universe model with isotropic perturbations 
up to the linear order. Hereafter, we call this model the inhomogeneous and isotropic 
universe model. By adopting the Newtonian gauge (see e.g., 20]), the infinitesimal world 
interval is written in the form, 


ds 2 = - [1 + 2^(t, x)] dt 2 + a 2 (t) [1 - 2 ip(t, x)] [d\ 2 + S 2 (x 5 k)dQ 2 ] , (2.1) 
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where a(t) is the scale factor scaled so as to be unity at present time t = t 0 , dil 2 is the 
line element of the unit 2-sphere, and, denoting the spatial Ricci curvature scalar of the 
background at present time by 6k, the function S(x j k) is defined as 


s(x;k) 



( 2 . 2 ) 


We assume that this universe model is filled with non-relativistic matter and the cosmological 
constant A, i.e., the so-called ACDM model. The stress-energy tensor of the non-relativistic 
matter, i.e., dust is given by 


Tau = p(l + 5)uM v + P (uju„ + U^Uu) , 


(2,3) 


where p and u M are the energy density and the 4-velocity of the background, respectively, 
whereas 5 and 5u M are the density contrast and the perturbation of the 4-velocity, respec¬ 
tively. The coordinate system is chosen so that the components of the background velocity 
4-vector and its perturbation are given by u ^ = (—1, 0, 0, 0) and Su M = (—ip, av x , 0, 0). 

The Einstein equations lead to the Friedmann equation for the background; 

2 


H 2 (a) := 


1 da 
a dt 


87rp 0 k A 


(2.4) 


3a 3 a 2 ' 3 ’ 

where p 0 is the background energy density at present time t = t 0 . Denoting the present 
value of H by Hq, the above equation is rewritten in the form 


h-‘ = Hi (’% + % + n A 

a 6 a z 


where 




8np 0 

3 HI' 


flk = ~ 


H 2 

n o 


and Da = 


A 

W 


The Einstein equations lead to the equations for the linear perturbations; 


k 


ip + ±HiP + [H z + 2- - — \ip = 0, 


(2.5) 


( 2 . 6 ) 


(2.7) 


5 


47rp 0 


^2 


d x ( S 2 d x ip ) + 3 kip 


-3/, 


( 2 . 8 ) 


= 


4vrp 0 


d x ip + Hd x ip 


(2.9) 


4 



where a dot denotes a partial differentiation with respect to t, and / is the velocity potential 
which is related to the perturbation of 4-velocity through 


v x = 


d xf 

aH' 


( 2 . 10 ) 


The general solution of Eq. (12. 7ft is represented by the linear superposition of the growing 
factor D + (t ) and the decaying factor D„(f), which are defined as 

aild < 2 - n > 

Hereafter, we assume that the decaying mode vanishes, since this assumption is consistent 
with the inflationary universe scenario. Thus, we have 


1 l’(t,x) = Hx)D+{t), (2-12) 

where h(y) is an arbitrary function of the radial coordinate y. Through Eqs. (12. 8p and (I2.9f) . 
5 and v x are expressed by using -0- Note that there are three parameters po, k, A and one 
arbitrary function h(y) i n the inhomogeneous and isotropic universe model. 


III. NULL GEODESICS IN INHOMOGENEOUS AND ISOTROPIC UNIVERSE 
MODEL 

Hereafter we drop higher-order terms in calculations without notifying that. All equations 
below are valid in the linear order of the perturbations. 

We assume that the observer in the inhomogeneous and isotropic universe model stays 
at the symmetry center y = 0 so that the observer recognizes the universe to be isotropic. 
In order to get the distance-redshift relation for the observer, we consider a past-directed 
radial null geodesic which emanates from the observer. By virtue of the isotropy in terms 
of the observer, both k e and k& should vanish. One of the non-trivial components of the 
geodesic equations is given by 

4- [a 2 (l - 2 hD+)k x \ + 2 a 2 (D + d x h) ( k x f = 0, (3.1) 

a\ 

where A is the affine parameter. Equation (13.11) determines k x , whereas the null condition 
determines k* in the manner 

k f = —a(l - 2 hD + )k x . (3.2) 
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Then we have equations for t and x as 


(3.3) 

(3.4) 


dt 

dX 


= k\ 


dx 7 y 
— = k x . 
dX 


The redshift z for the observer is given by 

(^lV)ls 


l + z = 


= —a (1 — hD + + v x ) k x , 


(3.5) 


(k^u^) | 0 

where subscripts s and o mean the quantities evaluated at the source and the observer, 
respectively, and we have chosen the affine parameter so that —[k^u^) 0 is unity. 

We rewrite the equations for the radial null geodesic in the forms appropriate for later 
analyses. Equation (13. ip is rewritten in the form, 


1 ^ = -2 (h - h dDj 


dt 

dz 


k x dz \ dt 

By differentiating the logarithm of each side of Eq. (13.5p with respect to z, we obtain 
1 dk x 1 


From Eqs. 


k x dz l + z 
and (13.7lh we have 
1 


tt _ h ^£±\ ^ , n dh dv x 
dt ) dz + dz dz 


H-h 


dDi 


1 + z ' V dt 
From the null condition (13.2p . we have 

dt 


dt dh dv y 

Tz +D+ Tz~ ~dJ = ' 


dx 


— + a (1 - 2 hD+) —j— = 0. 

dz dz 

We express the radial null geodesic as a function of z\ 


t = t{z) + 5t(z), 

X = x(z) + Sx(z), 


(3.6) 


(3.7) 


(3,8) 


(3.9) 


(3.10) 

(3.11) 


where the quantities with a horizontal bar represent the background solution. From 
Eqs. (13.8p and (13.9h . we see that the background solutions t(z) and x( z ) satisfy 

1 


dt 

dz 


(1 + z)H' 


dx 

dz 


1 

H' 


(3.12) 

(3.13) 
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Here note that 


H — Ho\jQ, m (\ + z) 3 + Hk(l + z) 2 + Ha • 


(3.14) 


Equation (j3.8[) leads to the equation for the linear perturbations as 


H 


cl5t dH dD + , ^ dh dv Y 

+ -j-St - —^h + D + — - = 0, 

dz dz dz dz dz 


(3.15) 


whereas Eq. (13.9p leads to 

dSt 1 / dSx r 2 D+h \ 
dz 1 + z \ dz H ) 


0, 


(3.16) 


where we have used Eqs. (j3. 12[) and (13.131) . If the inhomogeneous and isotropic universe 
model is completely fixed, Eqs. (13.15(1 and (13.161) determine 5t and Sx- 

The boundary conditions for Eqs. (13.121) . (13.131) . (13.151) and (13.161) are given by 


t(0) = to and y(0) = 5t(0) = hy(0) = 0. 


(3.17) 


Note that h, D + and v x are related to each other through Eq. (j2.9[) in the manner 

2 H 2 


v x~ 


+ z) 


, . dD+ 


f=0. 

dz 


(3.18) 


IV. THE EFFECTS OF INHOMOGENEITIES ON THE ESTIMATE OF DARK 
ENERGY 


In cosmology, the observable domain is restricted mainly to the null hypersurface gen¬ 
erated by past directed null geodesics emanated from the observer. This restriction causes 
the difficulty in recognizing isotropic inhomogeneities, since it is extremely difficult without 
changing the observation site over the cosmological scale to check whether such inhomo¬ 
geneities exist. As mentioned, we study the possibility that the effects of these formidable 
isotropic inhomogeneities and the functional degree of freedom of w in the equation of state 
of the dark energy (11.11) degenerate. 

We assume that the hypothetical observational data determines the dependence of the 
angular diameter distance d\ on the redshift z in the form 

dA = D a (z). (4.1) 
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Hereafter, we focus on the situation in which if we assume general relativity and the FLRW 
universe model, the observational data (14. ip indicates that the universe is dominated by the 
non-relativistic matter and the dark energy whose w in the equation of state (11.11) is written 
in the form 


w = ^2w n (l - a) n . (4.2) 

n =0 


where a is the scale factor, and w n is constant and w 0 < —1/3. We assume that w n = 0 
for n > 2. Furthermore, for simplicity, we assume that the observationally indicated FLRW 
universe model has flat space, k = 0. This FLRW model is characterized by four parameters, 
Hubble constant H 0 ,w 0 , w± and the present value of the energy density pdo °f the dark energy, 
or equivalently, 


fin : — 


87rp d o 

3 Hi ' 


(4.3) 


In the case of the real observation, the parameters H 0 , Wo, W\ and fid are determined so 
that the observed distance-redshift relation (14.11) is well fitted. However, the present analysis 
is theoretical, and hence instead of fixing the functional form of hypothetical observational 
data D A (z), we determine the cosmological parameters H 0 , w 0 , uq and fid first. Then we 
obtain the angular diameter distance of this FLRW universe model and identify it with 
D k {z). 

As mentioned in Sec. HI our purpose is to study whether the hypothetical observational 
data (14.ip can also be explained by the inhomogeneous and isotropic universe model given 
in Sec. m within the framework of the linear perturbation theory. In this case, the basic 
equations for the radial null geodesic (13.151) and (I3.16P and (I3.18P should be regarded as 
the system of differential equations to determine the inhomogeneous and isotropic universe 
model, i.e., the parameters Ho, fl m , flk, flv and the functional degree of freedom h(y); We 
will rewrite them into the forms appropriate for this purpose. 

In the inhomogeneous and isotropic universe model, the angular diameter distance from 
some light source to the central observer is equal to the areal radius at which the light is 
emitted; 


d A (z) = a(t ) 1 - D+(t)h(x) S(x',k ) 


x=x(z) 


= a(t)S(x ; k) + a(t) [(HSt - D + h)S(x', k ) + S"(y; k)Sx] , (4.4) 
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where a prime means a derivative with respect to its argument. Then, the angular diam¬ 
eter distance (14.41) is assumed to satisfy Eq. (14.11) . The angular diameter distance in the 
background FLRW universe model d A (z) is given by 


We define 


d A (z) 


s(x;k) 

1 + z ' 


Sd A (z) = D a (z) - d A (z). 


(4.5) 


(4.6) 


Then, Eq. (j4.4[) leads to 


8x{z) 


(1 + z)5d A (z) - \H(z)5t(z) - D+(t)h(x )] S{x\ k) 
S'(x;k ) 


(4.7) 


Equation (14.71) is equivalent to the hypothetical observed distance-redshift relation (14.11) . 

We regard Eq. (13.18f) as the differential equation for h. By substituting Eq. (14.7ft into 
Eq. (13.16ft . we eliminate Sx from Eq. (13. 16j) and obtain the differential equation for 5t. In 
order to eliminate dSt/dz and dh/dz from Eq. (j3. 15jl . we use Eq. (j3. 16[) with 8\ eliminated 
and Eq. (13.181) . and then, we obtain the differential equation for v x . As a result, we obtain 
the following system of differential equations to determine the inhomogeneous and isotropic 
universe model from the hypothetical observed distance-redshift relation (14.11) : 

^ = H(v x ,z), (4.8) 

d5t AT (h, 5t,v x , z) 

dz 

dVy , , . , r 

— = V (ft, St, v x , z ), 


(4.9) 

(4.10) 


where 


n 

V(z) 


J\f (h, St , v x , z) 


3i7gfl m (l + z) 2 v x 
2 H 2 


D + — (1 + z) 


dD 4 
dz 


(4.11) 


1 - 


H S(x; k) 
(! + «) S’fakY 


(4.12) 


1 + z 


d f {1 + z)Sd A 


dz 


S' 


+ 


- 1 + d( s JL' ]Ut 


dz V S' 
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(4.13) 


2D A 


H 


d 

dz 


D+S 

~S r 


h 


D+S 


S' 


n 


v = 


dD -1 
dz 


■h 


D + n + ^-6t + H^-. 

dz V 


(4.14) 


The boundary conditions for Eqs. (I4.8p - fl4.10p are given as follows. The value of h\ z=0 
can be made zero by the rescaling of the coordinate and introduce a new scale factor, (see 
Appendix A). Hence we impose 


h\z=o = 0. (4.15) 

The boundary condition on 5t(z) is given by Eq. (13.171) . Imposing C 1 regularity for h at the 
origin, that is, d x h |~=o = 0, from Eq. (12.911 . we obtain 

v x \ z =o = 0. (4.16) 

We numerically solve Eqs. (I4.8I1 - (I4.10I1 under the boundary conditions (j3. 17f) . (14.1511 and 
(14.1611 . If we find numerical solutions are so small that the linear approximation is valid, 
we may say that the FLRW universe model filled with the dark energy whose equation of 

state is given by Eq. (14.2p with a set of w n can also be explained by the inhomogeneities 

and isotropic universe model with the dark energy of wq = —1 and w n — 0 of n > 1, i.e., 
the cosmological constant. 


V. NUMERICAL PROCEDURE 

Before numerically integrating Eqs. (I4.8[) — (j4. lOj) . we need to fix three cosmological pa¬ 
rameters Ho, Ok and Ha of the inhomogeneous and isotropic universe model. The Hubble 
parameter H 0 should be the same as that of the FLRW universe model giving Da{z), since 
it is an observable quantity. By contrast, as explained below, one of f2 k and Ha cannot be 
freely determined as long as we demand that the solutions are at least everywhere C 1 . 

Equation (14.91) shows that this equation has a regular singular point at z — z cr which is 
a root of V(z) = 0. The function A f(z) := J\f(h(z), 5t(z),v x (z), z) should satisfy 

A f(z cr ) = 0 (5.1) 

so that the solutions of Eqs. (j4.8jl — (j4. 10H have finite and continuous derivatives with respect 
to z at z = z CT . The condition (15. ip leads to a relation between H k and Ha- In this paper, we 
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fix f^A = 0.7, and then we search for which guarantees the smoothness of the solutions 
at z — z cr by a kind of the shooting method as shown below. 

Since we can find z CT > 1 for the cases of our interest, we solve Eqs. (j4.8j) — (14. 10p from z = 0 
to z = 1 by imposing the boundary conditions (13.171) . (j4. 1 5j) and (14.1 bp and, at the same 
time, from z = z cr to z = 1 , by making a guess at h\ z=Zci and v x \ z=Zci and then fixing 
St(z cv ) so that Eq. (15.ip is satisfied. If we fail to get smooth solutions of Eqs. (I4.8I) - (I4.10I) . 
we select different values of f2 k , h\ z=Zci and v x \ z=Zci in accordance with the Newton method 
in the three-dimensional parameter space and then again integrate Eqs. (j4.8p — (j4. 10p from 
z = 0 to z = 1 and, at the same time, from z — z CT to z — 1. We iterate this procedure until 
we obtain numerical solutions sufficiently smooth at z — 1. Next, we integrate Eqs. (14.8H 
(14.10p outward from 0 = z CT with the values of f2 k , h\ z=Zci and v x \ z=Zci which guarantee the 
smoothness of the solutions. 

VI. NUMERICAL RESULTS 

Here, we introduce the density contrast in the synchronous comoving (SC) gauge A 
defined as 


A := 5 + 3/. 


(6.1) 


where / is the velocity potential introduced in Eqs (12. 8p and (12. 10p . The gauge transforma¬ 
tion from the Newtonian gauge to the SC gauge is given in Appendix B. 

At present, we do not have any observational data of the distance-redshift relation in 


the domain of z > 2, except for positions o 
Microwave Background Radiation (CMBR) 


le acoustic peaks in the spectrum of Cosmic 


21|. Hence we assume that the universe model 


obtained by solving the inverse problem is available in the domain of z < 2 only, and the 
domain of z > z^ {z\, > 2) agrees with the FLRW universe model with k — 0, i.e., the 
concordance ACDM model, since the concordance ACDM model is consistent with many of 
inflationary universe models. We call this concordance ACDM universe model the faraway 
background. It should be noted that the faraway background does not necessarily agree 
with the background universe model introduced in the previous section, and hence we have 
to introduce the density contrast relative to the faraway background. 
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We assume that the present energy density satisfies 

rrb 

/ A f (r 0 , r)r 2 dr = 0, (6.2) 

Jo 

where Af(ro,r) is the density contrast relative to the faraway background at the present SC 
time r = To, and r b is the comoving radial coordinate at which z = Z\> (see Appendix C on 
the SC coordinates r and r). The condition (16.2ft is necessary so that the domain of r > r b 
is exactly the concordance ACDM model. But we do not impose any additional conditions 
on A f of 2 < z < besides Eq. (I6.2|h 

In Figs. HM we depict the present density contrast A f | T=T0 in the homogeneous and 
isotropic background with fA = 0.7 as a function of the redshifht z, which can explain the 
angular diameter distance agrees with D\(z) of the FLRW universe model with the dark 
energy of various hid, u>o and W\. The outermost data points correspond to that of z = 2. 



FIG. 1: We depict the present gauge invariant density contrast Af| T=TQ in the homogeneous and 
isotropic background with Ra = 0.7 as a function of the redshift z, which can explain the distance- 
redshift relation of the FLRW universe model with various fid of the dark energy with wq = — 1 
and w\ = 0. 


In Fig. |U we depict Af| r=ro in the case of wq = — 1, w\ = 0 and various fid- In Fig. [2] 
we depict the same as Fig. [Tjbut in the case of fi d = 0.7 and Wi = 0 with various U7 0 - In 
Fig. [3j we depict the same as Fig. |T] but in the case of fi d = 0.7 and w 0 = —1 with various 
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FIG. 2: The same as Fig. [T] but fid = 0.7 and w\ = 0 for various wq. 



FIG. 3: The same as Fig. |1] but fid = 0.7 and wq = — 1 for various w\. 

w\. We can see form these figures that the isotropic inhomogeneities may mislead us about 
the estimate of the values of Wo and w±. 
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VII. PROBABILITY OF OCCURRENCE 


We have seen that the isotropic perturbations can cause systematic errors in the estimates 
of the cosmological parameters and the equation of state of the dark energy. 

In this section, we show the probability of occurrence of these isotropic perturbations. 
We assume that the initial density perturbations obey random Gaussian statistics which is 
consistent to the inflationary universe scenario: the Fourier transform of the density contrast 
is given by 

A(t, k) — j d 3 x A(r, x)e~ lk ' x , (7.1) 

and then, we have 

(A (t, k)A*(t, fc')> = (2 tt) 3 <^(fc - k')P(t, k), (7.2) 

where (• • •) represents the ensemble mean, S-£>(x) is Dirac’s delta function, and P(t, k) is the 
power spectrum at time t. We assume 

P(t, k) = A 0 k n T 2 (k)D 2 + (t ) (7.3) 


where n and A 0 are constants determined by observing the fluctuation of the CMBR, whereas 
T(k) is the matter transfer function; We adopt n and A 0 derived from the Planck data and 
the matter transfer function derived by Eisenstein and Hu 22]. The matter transfer function 
is determined by background cosmological parameters. 

The mass contrast within the radius r is given by 


SM, . 1 

- r, r = — 

M y V 


d 3 x d (r — |tc|) A (t, x), 


where 'd(z) is Heaviside’s step function, and 


(7.4) 


V — d 3 x ‘d (r — \x\) = -7rr 3 . 


(7.5) 


The mean square of the mass contrast is then given by 

2' 


a r, r := 


SM 

~M 


(r,r) 


27T 2 


P(t, k)W(kr)k 2 dk, 


'o 


where 


9 


W(y) = —(sinj/ - ycosj/)' 


(7.6) 


(7.7) 
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FIG. 4: IV max is depicted as a function of II c j We fix the parameters of the FLRW universe with 
dark energy as Ha = 0.7, wq = —1.0 and w\ = 0. 
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FIG. 5: The same as Fig. SI but = 0.7 and w\ = 0 for 


various Au>o- 


The present mass contrast of the inhomogeneous and isotropic universe model obtained 
by solving the inverse problem is defined as 


A m W : = 


Y f A { (T 0 ,y)y 2 dy 


(7.8) 
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iV, 



FIG. 6: The same as Fig. 01 but Sid = 0.7 and wq = —1 for 


various w \. 


Then, we define 


N{r) := 


A m (V) 


(7.9) 


ff(To,r) 

The larger N{r ) implies that such an inhomogeneity is rare. The maximal value of N{r ) in 
0 < r < r 2 is denoted by fV max , where r = r 2 is the comoving radial coordinate of £ = 2. 

In Figs. HHH we depict the fV max with the dark energy of various fid, Wq and w\. In 
Fig. 01 we depict fV max as a function of fid in the case of wq = —1, w\ = 0. fid close to 
0.7 means lower density contrast. In Fig. 01 we depict the same as Fig. 0]but in the case of 
fid = 0.7 and w\ = 0 with various wq. In Figj6l we depict the same as Fig. 0]but in the case 
of fid = 0.7 and w 0 = — 1 with various w\. 


VIII. SUMMARY AND DISCUSSION 

We studied the systematic error caused by isotropic inhomogeneities, which appear in 
determining the amount of the dark energy and its equation of state. We have provided 
the hypothetical observational data which are equivalent to those obtained in the universe 
dominated by the dark energy of various fid and equation of state. Then we solved the 
inverse problem to evaluate the density contrast in the FLRW background filled with non- 
relativistic matter and the cosmological constant, so that the hypothetical observational 
data is explained. The probability of occurrence of these isotropic perturbations is also 
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evaluated in accordance with the standard inflationary universe scenario together with the 
Planck data. 

Our result implies that the uncertainty in Wq is ±4 x 10~ 3 in 10cr confidence level, and even 
if we get the estimate w 0 = —1.005, it does not necessarily mean the existence of the phantom 
energy. Even in the case that observational data suggest w 0 = —1.01, it is impossible to 
deny the possibility that it comes from a large scale isotropic inhomogeneity; since it may be 
only one realization within our observable domain, we can not exclude, from the statistical 
argument, the possibility that such an inhomogeneity unexpectedly appears. Although it 
is very challenging to observationally determine how large the isotropic perturbations are, 
we need constraints on their existence in order to arrive at the conclusion with sufficiently 
confidence. This is the future problem. 
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Appendix A: Rescaling to eliminate the metric perturbation at the symmetry center 

In this paper, we consider the isotropically perturbed FLRW universe model filled with 
non-relativistic matter and the cosmological constant. As shown in Sec. Qdl the metric is 
given by 

ds 2 = - [1 + 2 h(x)D + (t)} dt 2 + a 2 (t) [1 - 2 h(x)D + (t)] [d\ 2 + S 2 (x; k)dfl 2 ] . (Al) 

We rescale the time and radial coordinates in the manner 

dt = [1 + h(0)D + (t)]dt, (A2) 

X = [1 - h(0)D + (t o )]x, (A3) 

and introduce a new scale factor defined as 

d(i) = [1 - h( 0) {D + (t) - D+(t 0 )}] a(t). (A4) 
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Then, up to the first order perturbations, the line element is rewritten in the form 

ds 2 = - [l + 2 h(x)D + (t)\ dt 2 + d 2 (t) [l - 2 h(x)D+(tj\ [dx 2 + S 2 (x; k)df2 2 ] , (A5) 

where 


k = [1 + 2h(0)D + (t o )] k, (A6) 

Hx) = h(x) - h{0), (A7) 

D + (i) = D+(t). (A8) 


Note that h vanishes at the symmetry center y = 0 = X- 
The Hubble equation for the new scale factor a is given by 


1 da(t)\ 


a(t) dt 


I = [l-2h(0)D + (t)} 


1 da(t)\ 


) 


= 


= Hi 


_a 3 {t) 


a(t) dt 
l + 3h(0)D + (t) - 


- 2/i(0) 

2h(0) \ 


1 da(t)\ dD + (t ) 


a(t) dt J dt 


n m H 2 ) 


Hi< 


+ [1 + 2h(0)D+(t)] 9 . + Ha 


a 2 (t ) 


+ 




+ Ha 


(A9) 


0 \d 3 (t) ' a 2 (t) 

where we used the original background Hubble equation (12.5ft and Eq. (12. lip and furthermore 
introduced the following new cosmological parameters: 


fi m 

Hk 


H m 


1 + 3h(0)D + (to) — 


Hk [1 + 2h(0)D + (to )] • 


2h(0) 
H W 2 


(A10) 

(All) 


Equation (1A9K implies that a is regarded as the scale factor of the FLRW universe model of 
the cosmological parameters H m , Hk and Ha- Hence, without loss of generality, we can get 
the isotropic metric perturbations to vanish at the symmetry center x = 0 by choosing an 
appropriate background FLRW universe. 


Appendix B: Gauge transformation 


The transformation between the Newton and synchronous comoving (SC) gauges is given 


by 


dr = (1 + ijj)dt — av x dx, (Bl) 

dx = ——dt + f 1 + a _ d Y w ) dx- (B2) 

a V 4vr Po J 
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By integrating the above equations, we have 

t = t+ {j’ ( B3 ) 

* = X + (B4) 

where we have used Eq. (12.9p in obtaining r, whereas Eq. (12. 10p has been used in obtaining 
r. The infinitesimal world interval in the SC gauge is given by 

ds 2 = — dr 2 + a 2 (r) [1 — 2(r/> + /)] 


x 


) dx 2 + S 2 (x] k) (1 — ' j} dx'ijj ) di ? 2 


2ttpo 1 / V 2vrpo S(r;k) 

By the transformation (IBID and (1B2|) . we have, up to the first order, 


dt\ 


dr , 


dt \ f dx 


dr J ^ \dr 


Xr = (^) T it + 2(^] (-£) T tx +(-£) T«=p(r)(l + i + 3/), 


d\ 


where T tt = p(t)(l + 5 + 2'ijj) has been used. Then, since T tt = p(r)(l + A), we have Eq. 


(B5) 


(B 6 ) 

(ED- 


Appendix C: Replacement of the background 


Here, we consider the replacement of the background from the FLRW universe model of 
k 7 ^ 0 to that of k — 0. The FLRW universe model is assumed to be filled with the dust 
and the cosmological constant. 

We write the FLRW universe model with k ^ 0 in the form 

ds 2 = —dr 2 + a 2 (r) ^1 + ^R 2 ^ (dR 2 + R 2 dVt 2 ) . (Cl) 

The relation between the radial coordinates x dr Eq. (12. ip and R is given by 


R = 


1 — \J 1 — kS 2 (x ; k ) 


(C2) 


kS(x;k ) 

We assume \kR 2 \ < \k/po\ -C 1. Then, Eq. (ICip is rewritten in the form of the infinitesi¬ 
mal world interval of the FLRW of k = 0 with perturbations as 

k 


ds 2 = -dr z + A 2 (r) 


1 + 2A a - -R 2 

a 2 


(dR 2 + R 2 dQ 2 ) , 


(C3) 


where A(r) is the scale factor of the FLRW universe model with k = 0 and 

a(r) - A(r) 


A„ := 


A(r) 


(C4) 
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with the assumption |A a | -C 1. By the definition, we have 


1 da 
a dr 

On the other hand, from Eq. (1C4|) . we have 


87rp 0 k A 
3a 3 a 2 3 


a dr J 

We may assume that A satisfies 


1 da\ 2 /I dA\ 2 2 dA dA a 


A dr J A dr dr 


1 dA\ 87Tp 0 A 
A dr J 3A 3 3 

Then, we have, from Eqs. (1C5|) and flQ6[i . 

8np 0 A k 2 dA dA a 
A 3 a+ A 2 + A dr dr ~ ' 

Here it should be noted that Eq. (12. 11 j) leads to 

2 k\ . , 2 - d(aD + ) 

8np -- ) (aD + ) - - + 2 H -—— 


a“j or dt 

By comparing Eq. (1C8I) with Eq. (lC9jl . we find 

A a = 


= 0 . 


(C5) 

(C6) 

(C7) 

(C8) 

(C9) 

(CIO) 


where 

np 4 ) p dB 

+ A J B 3 H 3 (B) 

with 

»» ; l/WI 

The density contrast of the FLRW with k ^ 0 relative to that of k = 0 is defined as 

A FL rw := (§•)“' (| - f) = -3A„ = \kA(r)V + (T). (Cll) 

Note that A is not unity at present time r = To, if a is. By the definition, we have 


A ( r o) — 1 — A a0 j 

where A a0 := A a (ro). We introduce the following quantities; 
a = (1 + A a0 ) A, r = (1 — A a0 ) i?, k — (1 + 2A a0 ) k and po = (1 + 3A a0 ) po- (C12) 
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Then, the infinitesimal world interval becomes 


ds 1 2 = —dr 2 + a 2 (r) [ 1 + 2A a — ^-r 2 ) (dr 2 + r 2 dfl 2 ) , 


(C13) 


The density contrast A F lrw is written as 


Aflrw = ^kA(r)V + (r) = ^kd(r)D + (r), 


(C14) 


where 


with 


_ H(d) r dB 
a J B 3 H 3 (B) 


H(B) 




By the result obtained here, we can derive the density contrast Af relative to the FLRW 
universe model with k — 0 from the density contrast A relative to that with k ^ 0 in the 
manner 

Af = A + Aflrw = A + -fca(r)l) + (r). (C15) 

By using the result obtained in this appendix, the infinitesimal world interval (IB HI) may 
be rewritten in the form 


where 


ds 2 = 


-dr 2 + d 2 [ 


t 


+ /) 


x 



dr 2 + r 2 



a i q ? 
--Orlp 


2vrp 0 r 



(C16) 


7 , . 27tpo^ 2 Aflrw 

^ = il> +-—-, 

Sa 

r r 27rp 0 r 2 A F LRw r 2 Aflrw 

3a 6 dD + 


(C17) 
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